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ABSTRACT 

Non-commutative geometry at inflation can give arise to parity violating modulations of the pri- 
mordial power spectrum. We develop the statistical tools needed for investigating whether these 
modulations are evident in the Cosmic Microwave Background (CMB). The free parameters of the 
models are two directional parameters (9,4>), the signal amplitude A.*, and a tilt parameter n* that 
modulates correlation power on different scales. The signature of the model corresponds to a kind of 
hemispherical power asymmetry. When analyzing the 7- year WMAP data we find a weak signature 
for a preferred direction in the Q-, V-, and W bands with direction (I, b) = (-225°, -25°) ± (20°, 20°), 
which is close to another previously discovered hemispherical power asymmetry. Although these re- 
sults are intriguing, the significance of the detection in the W-, V- and Q-bands are nonzero at about 
2ct, suggesting that the simplest parameterization of the leading correction represents only partially 
the effects of the space-time non-commutativity possibly responsible for the hemispherical asymme- 
try. Our constraints on the presence of a dipole are independent of its physical origin and prefer a 
blue-tilted spectral index n* ~ with the amplitude A.* ~ 0.18. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

During recent years, studies of the cosmic microwave 
background (CMB) have greatly improved our under- 
standing of the early universe. Observations of the 
CMB anisotropics, such as those obtained by the Wilkin- 
son Microwave Aniso t ropy Probe (WMAP ) experiment 
([Bennett et all 120031 : iHinshaw etall 120071 ). have pro- 
vided us with a deep insight on the composition of 
structure and energy in our universe, giving rise to 
the ACDM model. The ACDM model requires that 
the universe has undergone an epoch of rapid acceler- 
ated expansion. This epoch is named inflation, and 
is thought to h ave been driven by a single or sev- 
eral scalar fields (jGuth et all 119811 ). In addition, the 
model of inflation establishes a highly successful theory 
for the formation of primordial density perturbations, 
providing the required seeds for the large-scale struc- 
tures (LSS). These large scales structures are observed 
in the CMB today (iGuth et al l 119811: iLinde et alJ | 1982 



Mukhanov fc Chibisovl 119811 iStarobinskv et al I Il982 



Linde et all I1983L 119941: iSmoot et al 1119921 : iRuhl et al 



20031: iRvnvan et all 120031: iScott et al.H2003l) 

Inflation can explain why the observed universe 
should be nearly isotropic on large scales. How- 
ever, anom alies found in the CMB during the re- 
cent y ears (Ide Oliveira-Costa et al.l 120041: iVielya et all 



200l"lEriksen et all l2004at iGroeneboom et all 120101 : 
Hoftuft et al.ll2009[ ) suggest that some anisotropy could 
be present at inflation too. 

Several theoretical possibilities were put for- 
ward t o explain su ch ano ma lies (| Ackerman et all 
2001 iKanno et al. I [2001 iDimopoulos I 



Yokovama and Soda I I2008T Barrow and Hervik I 12009 : 
Esposito-Farese.Pitrou and Uzan I 120091) . One of the 



rotational invariance in the early universe. The presence 
of such vector fields lead to quadrupole modulations 
of the CMB anisotropy, which fo r vector perturbations 
was shown b vlDurrer et all (119981 ) . Rec ently, it was also 
suggested by IGroeneboom et all ([20101 ) that the 5-year 
WMAP data contains a significant signal of a primordial 
vector field which would break rotational invariance, 
corresponding to a 9er detection in the W-band . How- 
ever, several authors have since then claimed that the 
signal is due to systematic effects, as the rotational axis 
is aligned with the rotational axis of the satellite. A 
likely candidate for the source of t he signal is th e refore 
asymmetric instrumental beams (Hanson et al.l 120101 : 
Komatsu et~all 120101: iBennett et all 120101 ) . Inaddition. 



Pullen fc Hiratal (j2010f l have shown that there is no sig- 



nal of breaking of rotational-invariance in the observed 
LSS data, which makes less likely that the signal is of a 
cosmological origin. 

Unfortunately, su ch scenarios seem to be plagued 
with ghosts ([Himmetoglu. Contaldi and Peloso I 
l2009a| ) and several other sorts of insta- 
bilities dKoivisto. Mota and Pitrou I 



Him metoglu. contaldi and Peloso 1 
Himmetoglu. contaldi. Peloso I l2009bl ). 



but 




generalizations and alterna t ives can also be considere d 
Boehmer and Mota I (|2008l) ; IKoivisto and Mota] (|2008j) : 



possibilities is to introduce a vector field that breaks the 



Jimenez et al. I (|2009l i: lArmendarizPicon I (|2007f ? . 

Several of the anomalies, such as the axis of evil or the 
hemispherical asymmetry in the power spectrum, raise 
the question whether there could exist odd-parity mod- 
ulatio ns in the CMB. It was shown by IKoivisto fc Motal 
(|2010f) that the effects of space-time non-commutativity 
at inflation can generate modulations. The leading con- 
tribution, as generically in any parity-violating model, is 
a dipole modulation, but up to now it hasn't been explic- 



2 



itly looked for from the observed microwave sky. There 
is a general framework for describing the effects on the 
CMB anisotropy in terms of bipolar harmonic expansion 
(lHaiian and Souradeep I 120031 : iPullen fc Kamionkowskil 
I2007H which we we apply here. The theoretical back- 
ground is discussed in section [21 while the numerical im- 
plementation is discussed in section [3] In section 01 we 
investigate the model properties and conclude that the 
effect is similar to a hemispherical dipole asymmetry. We 
perform several analyses on simulated data and discuss 
the joint posterior results. In section [5l we analyze the 
seven-year WMAP sky maps and find weak evidence for 
a preferred direction located at the previously describe d 
hemispherical power asymmetry bv lHoftuft et al.l (|2009f ). 
In [6] we conclude our results, and discuss what how more 
general models might give rise to more significant results. 

2. BACKGROUND 

The parity violating effect on the primordial power 
spectrum studied in the present paper can be generated 
by quantum effects during or shortly after inflation. Such 
effects have yet to be considered due to the popular as- 
sumption that the power spectrum should be invariant 
under spatial inversion. In principle there are no phys- 
ical arguments that exclude them as long as the physi- 
cal observables generated from the theory are consistent, 
so one should therefore test this class of models against 
WMAP data and assess their viability. The odd-parity 
bispectrum was recently explored in a similar spirit by 
iKamionkowski and Souradeep I (120 lOD. 

It was shown in iKoivisto fc Motal (|2010h that by taking 



into account non-commutativity of space-time at very 
short distance scales, the inflationary power spectrum 
becomes direction-dependent. This result can be derived 
in several ways. The canonical 0-type non-commutation 
relations are 

(1) 



where 9^ v is an antisymmetric constant matrix of di- 
mension length squared that in the simplest case is a 
constant in a given coordinate system. Effectively this 
results in a non-commutative field theory. Another way 
to incorporate (stringy) uncertainty principle into geom- 
etry is to define a deformation of the Heisenberg algebra 
of the quantum fields. It was shown that with simple as- 
sumptions, in both approaches one may predict, a power 
spectrum of the forn{3 



P 9 (k) - P (k)e me k (2) 

where Po(k) is the (isotropic) power spectrum when 
the non-commutativity is turned off. H is the Hub- 
ble factor and 9 is the time-space component of the 
non-commutativity matrix, which sets the scale for the 
strength of the dipole modulation correction, or Ai m ~ 
\9\ in the present formalism. For a com parison with other 
literat ure see the references cited by IKoivisto Sz Motal 
(|2OT0h . The main feature of this model is the parity- 
violating power spectrum referenced in eq. © which 

1 More generally, one obtains Pg(k) = Po(k) [a cos(H8 ■ k) + 
i/3 s'm(H9 ■ k)] where a and /3 are parameters typically of order one, 
and the trigonometric functions can be replaced by hyperbolic ones 
depending on the prescription. However, at the leading order the 
behavior is qualitatively the same. 



implies that the two-point correlations, or the covariancc 
matrix, generates consistent physical observables such as 
CMB maps without any further modifications. 

However, it may turn out that similar criticisms apply 
also for these early universe models as were mentioned 
against models which br eak rotation invariance at early 
times such as the one bv lAckerman et all (|2007t ). In par- 
ticular, the present realisations of field theories involv- 
ing space-time non-commutativity may not be unitary. 
As a fully consistent underlying non-commutative the- 
ory is still lacking, obviously we cannot make a unique 
and sharp prediction of the fine details of its implica- 
tions to cosmology. They depend upon the inflationary 
model and the form of non-commutativity (for example, 
whether 0^ is constant in Eq.([T]) and in which coor- 
dinates). These objections notwithstanding, our search 
of a dipole signature is strongly motivated since it is 
a generic prediction of these theories, and it seems a 
quite unique prediction too, since it appears to be oth- 
erwise difficult to realize the situation where the primor- 
dial spectrum of fluctuations is not symmetric. It is 
worth stressing also that order-of-magnitude estimates 
show that the theoretically plausible energy scales of 
non-commutativity are already probed by the present 
CMB experiments and the ensuing bounds will be pushed 
higher by the PLANCK satellite data. Thus, testing the 
leading order signature (dipole), gives us clues what to 
expect in a given non-commutative model for the next- 
to-leading signatures (higher multipole modulations and 
non-gaussianity). This provides an observational handle 
for the construction of viable and consistent theory. 

The prediction for the CMB anisotropy pattern in non- 
commutative geometry begins with the expansion of the 
primordial spectrum: 
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where the para m eteriz ation by 

lArmendariz-Picon fc Pekowskvl (|2009t) is employed. 
The sum over L could in principle cover the entire 
spectrum of multipoles and M runs from — L to L. The 
resulting signal covariance is written as 
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where the spectral-index tilm that parameterises the 
scale-dependence is included in the integrated contribu- 
tion over all scales, 
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e e (k)e t ,(k). (5) 



The multipole moments of the sources 0/>( fc) are com- 
puted using a modified version of CAMB (jLewis et al.l 
I2000D with the spectral index as an input parameter. The 
geometrical facto r s tim -i'm 1 in eq. (@| arc provided by 
IKoivisto fc Motal IpOloT) . The non-commutative nature 
of the fields responsible for the perturbations give rise to 
non-hermitian signal covariance Sp ,. tm — —Sg m; £i m i for 
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the dipole (L = 1) contribution which is the only term 
in the expansion that is considered in the present paper. 
The naive expectation from both theory and observa- 
tions is that higher order terms are suppressed. This 
make this assumption here but it should be examined 
more thoroughly. The anisotropy in Q is then added 
to the isotropic signal, corresponding to the L = term. 
One obtains a direct interpretation of Aqq = A as the 
primordial scalar amplitude in the canonical expression 
for the isotropic matter spectrum. Furthermore 



Aim = Aim?m 



(6) 



where Aim is the amplitude which we are estimating. 
The constant of proportionality in all three cases is 
i/\/3rkoA where r = \tq\ = H9. Here \[9 is the non- 
commutative lcngthscale, fco is the pivot-scale. The unit 
vectors appearing in cq. ((6]) are the spherical vectors 
paramctcrising the direction of the anisotropy, 



r± = =F 



r x =F ir 
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(7) 



In the case |T]) with 9 constant in the comoving frame, we 
have 9 = \9\ and n\M = 2, but as argued above this need 
not be the case for general models. Thus we will check 
this particular case separately but otherwise keep both 
the spectral index and the amplitude as free parameters. 
So in all the model contains six unknown parameters, the 
three spectral indices uim, the non-commutative length- 
scale and a direction (9, 4>) contained in the unit 
vectors of cq. (0. However, typically the slope of the 
modulation does not depend upon the azimuthal orien- 
tation. Then the number of parameters is reduced to 
four: the direction (9,4>), a coupling A* and the reduced 
spectral index n* for the L = 1 case. 

3. METHODS 

We now discuss the method for mapping out 
the desired posterior. The method for obtaining 
the posterior is similar t o the method presented 
in IGroeneboom fc Eriksenl (|2009T l : iGroeneboom et all 
(|2009t ). with the exception of a new covariance matrix 
and a new parameter. CMB data observations can be 
modeled as: 

d = As + n, (8) 

where d represents the observed data, A denotes 
convolution by an instrumental beam, s(9,4>) = 
J2e m a tmYt m {9, (j)) is the CMB sky signal represented in 
either harmonic or real space and n is instrumental noise. 
It is generally a good approximation to assume both the 
CMB and noise to be zero mean Gaussian distributed 
variatcs, with covariance matrices S and N, respectively. 
In harmonic space, the signal covariance matrix is de- 
fined by Se mt ii m ' = (ae m cigr m ,) ■ In the isotropic case, 
this matrix is diagonal. The connection to cosmological 
parameters lo is made through this covariance matrix. 
Finally, for experiments such as WMAP, the noise is of- 
ten assumed uncorrelated between pixels, N^- = afSij, 
for pixels i and j, and noise RMS equals to cr,. 

Let oj denote a set of cosmological parameters. Our 
goal is to compute the full joint posterior P(w|d), which 
is given by P(u\d) oc P{d\u)P{u) = C(uj)P(uj) 7 where 



C(u}) is the likelihood and P(oj) a prior. For a Gaussian 
data model, the likelihood is expressed as: 



C(u>) (X 



e -id T c- 1 (^)d 



(9) 



where C = S + N is the total covariance matrix. 

3.1. The Gibbs sampler 

The problem of extracting the cosmological signal s 
and ui fro m the full signal b y Gibbs sampling was ad- 
dressed bvl Jewell et all (|2004l) . lWandelt et al.l(|2004H and 
lEriksen et al.1 (|2004bl ). The CMB Gibbs sampler is an 
exact Monte Carlo Markov chain (MCMC) method that 
assumes prior knowledge of the conditional distributions 
in order to gain knowledge of the full joint distribution. A 
significant fraction of the CMB data is completely domi- 
nated by galactic foreground, and about 20% of the data 
needs to be removed. This might sound trivial, but in 
reality it complicates processes as the spherical harmon- 
ics no longer are orthogonal. The Gibbs sampler solves 
this problem intrinsical ly, as the galaxy m ask becomes a 
part of the framework (|Groeneboomll2009t) . 

The main motivation for introducing the CMB Gibbs 
sampler is the drastically improvement in scaling. With 
conventional MCMC methods, one needs to sample the 
angular power spectrum, Cg = (ai m a^ m ), from the distri- 
bution P(Ce\d), which scales as C(iVp ix ), where N piK is 
the size of the covariance matrix. For a white noise case, 
the Gibbs sampler reduces this to O(N^). In other 
words, the Gibbs sampler enables effective sampling in 
the high-£ regime. 

3.2. Sampling scheme 

In order to sample from the full joint distribution 
P(Cti w,s|d) using the Gibbs sampler, we must know 
the exact conditional distributions P(s|CV,w,d) and 
P(Ce,uj\s). The Gibbs sampler then proceeds by alter- 
nating sampling from each of these distributions: 



-P(C t ,cj\**,d) 
-P( S \(C i ,w) i+1 ,d) 



(10) 
(11) 



The first conditional distribution is expressed as: 

-is T S( W )- 1 s 

P(C ilU \B,d) = —==-, (12) 

vpM| 

and is distributed according to an inverse Gamma func- 
tion with 11 + 1 degrees of freedom. The remaining con- 
ditional distribution is 



P(s\Ce, bJ, d) oc e 



-i(s-S) T (S(o;)- 1 +N- 1 )(s-S) 



(13) 



where s = N _1 d. In other words, P(s\C^,uj,d) is 
a Gaussian distribution with mean s and covariance 
(S(w) _1 + N -1 ) -1 . Numerical methods for sampling 
from t hese distributions were discussed bv iGroenebooml 

(pool . 

4. MODEL PROPERTIES AND PREDICTIONS 

In this section, we review the numerical setup of the 
analysis. Most of the framework is s imilar to the one em- 
ployed bv IGroeneboom fc Erikse n (2009). However, we 
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Fig. 1. — The non-commutative dipolc-corrclatcd covariance ma- 
trix S on logarithmic scales. Blue represents lower values (or zero), 
red contains larger values. 
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Fig. 2. — The transfer function Integrals 
of n*: 0.5 (black), 1.0 (red) and 1.5 (blue). The integral corre 
lates the power of the dipole signal for various scales I, and are 
normalized according to their area. Note that low values of n» 
corresponds to large-scale correlations and vice versa, large values 
of n* corresponds to small-scale dipole correlations. 

introduce a new parameter in addition to the anisotropic 
amplitude coupling A* oc Aoo and the directional param- 
eters 8 and <fi. The parameter n„ represents the spectral 
index of the correlation integral for the L = 1 case, and 
determines the tilt of the correlation power spectrum. 

4.1. Numerical setup 

We have previously developed a MCMC framework 
that enables sampling over sparse anisotropic universe 
models, meaning models predicting a covariance matrix 
that is sparse. This MCMC sampler is i ntegrated into 
a Gibbs samp ler named COMMANDER (jEriksen et al.l 
I2004U l2008al) . The Gibbs sampler alternates between 
sampling the CMB signal and the anisotropic model pa- 
rameters. The details o n how this method was imple - 
mented was described bv lGroeneboom fe Eriksenl (|2009f ) . 

The covariance matrix is for this model given as fol- 
lows: 



A-V L = °° L 



LAI 



L=l M=-L 



dipole modulations, or L = 1. A plot of the sparse co- 
variance matrix S can be seen in [T] 

The diagonal case L = represents the power spec- 
trum Ci, where we assume the tilt is equal to the spec- 
tral index ni m = n s from the best-fi t seven-year WMAP 
spectral index (jLarson et al.l l2010T l . We introduce the 
parameter n* to replace ni m for the L = 1 such: 



loot — 



dk 

~k 



e e (k)e t (k) (15) 



(14) 

Throughout this paper, we are only concerned with the 



where we have chosen fco = 0.05Mpc~ as the tilt scale. 
Three examples of I 1M for = 0.5, 1.0 and 1.5 are pre- 
sented in Figure [2] The integral in equation [15] needs to 
be pre-computed in numerical software such as CAMB 
(jLewis et al.ll2000D . It is difficult to implement CAMB 
into the existing Gibbs sampler framework, so we utilize 
a different scheme. First, we pre-compute 10 000 inte- 
grals for n* in an interval [—1,4] and store the data in 
a binary file. This interval is large enough to allow for 
all types of scale-dependence in the anisotropy, and is 
our prior for n*. For all purposes we consider ti* to be 
continuous. The precomputed binary file is then loaded 
into the anisotropic MCMC framework such that n* can 
be treated as a free parameter. 

When simulating CMB maps for this model, the con- 
nection to a hemispherical dipole asymmetry becomes 
imminent. This is depicted the top-left frame of Figure 
[31 where we present a simulated CMB map with A* = 3.0 
and n* = 1.0 with direction (l,b) = (224°, 22°). In or- 
der to verify that this really is a signal similar to a dipole 
modulation, we simulate a map with A* = 0.0 and divide 
the A* = 1.0 map with the A* = map. The resulting 
dipole structure can be seen in the remaining tree maps 
in Figure [31 where we have included simulations with tilt 
parameter n* = 0.5 and n* = 1.5. Note that from the 
integral in Figure [2] n* = 0.5 corresponds to a large- 
scale dipole modulation while = 1.5 contains more 
small-scale modulations. 

4.2. Analyzing simulated data 

We verify code by performing both a brute-force and 
MCMC analysis on a noise-free, unconvolved simulated 
CMB map with no sky cut. In order to build the joint 
two-dimensional distribution of n* and A*, we need an 
in creased amount of samp l es as compared to the case 
bv lGroeneboom fc Eriksenl ((2009). In addition, adding a 
new parameter will in general decrease the significance of 
the results if the parameters are correlated. As A* and n* 
can be expected to be correlated, we choose to normalize 
the integrals such that the area under each graph for all 
n* are the same. The degeneracy is broken, as can be 
seen in Figured) Here, we create a noiseless, unconvolved 
CMB map with ?i sidc = 128 and £ max = F^lf = 150 
adopting the best-fit ACDM model determined from 
the seven-year WMAP data (jLarson et alj[2010h . The 
model input parameters are = 1.0, A* = 1.5 and 
(l,b) = (224°, 22°). The joint and marginal posteriors 
are presented in Figure [U Note how the posterior of 
n* is distributed similarly to a x 2 -alikc distribution, and 
that A* and n* arc not degenerate. The distribution 
is always symmetric around A* = 0.0, due to negative 
amplitude A* in a direction (l,b) corresponds to a pos- 
itive amplitude in the opposite direction (— I, —b). This 
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Fig. 3. — The effect of the dipole-modulating model on simulated CMB maps. Top left: A simulated CMB map with a large dipole 
contribution. The remaning maps show simulated dipolc-modulated maps divided by a non-dipole-modulated map. Note how different 
values of the correlation tilt n„ induce large-scale or small-scale dipole correlations. 
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Fig. 4. — Results from a simple analysis of simulated data with- 
out foregrounds, noise or convolution. The input parameters were 
A* = 1.5, n, = 1.0 with direction (l,b) = (224°, - 22°). Top left: 
The directional posterior from a MCMC sampling. Top right: the 
directional posterior from a brute-force run. Lower left: The joint 
posterior P{A», n*|d). Note that these parameters are not degen- 
erate. Lower right: the marginal posteriors P(A*\d) and P(n* |d). 



is not a problem when A* is large when compared to its 
standard deviation, but for low values of A* it becomes 
difficult to separate the negative peak from the positive, 
as they will merge. 

We continue by creating realistic V-band differencing 
assembly (DA) simulations with n ns ide = 512, A* = 0.25, 
= 0.25 with direction (7,6) = (224°, -22°). The 
maps are produced using the V-band beam and noise 
properties. In addition, we add synchrotron, free-free 
and thermal dust foreground templates as described by 
iGold et al.l (|2010l ). The VI simulation is depicted in Fig- 
ure [5] The analysis is performed using the Gibbs sam- 



Simulated VI data, Gibbs sampler 
Directional posterior 



Simulated VI input data 
with foregrounds 




Fig. 5. — Results from a full-scale analysis of simulated realistic 
WMAP V-band data including foregrounds, noise and convolution. 
The input parameters were A* = 0.25, n* = 0.25 with direction 
(l,b) = (224°, -22°). Top left: The directional posterior P(n\d). 
Top right: The simulated VI map with foregrounds and a dipole 
modulation. Lower left: The joint posterior P(A*, n* |d). Note 
that these parameters are more degenerate for low values of .4*. 
Lower right: the marginal posteriors P(A*\d) and P(n„ |d). 

pier for ^"*° ff = 400 where we impose the KQ85 mask 
ijGold et all 12001) . which removes 18% of the sky. The 
analysis successfully reproduce the input parameters, as 
is seen in Figure [5] However, note that the distributions 
are wider than in the noiseless, perfect case, and that 
tail of the marginal posterior P(^4*|d) merges with the 
positive values of P(— A*|d) near zero. 

5. WMAP ANALYSIS 
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TABLE 1 

Summary of marginal posteriors from the seven-year 

WMAP DATA 



Band 


t 


range 


Mask 


Amplitude A* 


Direction (I, b) 


Wl-4 
Vl-2 
Vl-2 
Ql-2 


2 
2 
2 
2 


-400 

- 400 

- 400 

- 400 


KQ85 
KQ85 
KQ75 
KQ85 


0.17 ±0.08 
0.18 ±0.08 
0.18 ±0.08 
0.19 ±0.10 


(230°, 
(225°, 
(225°, 
(225°, 


-25°) ± 20° 
-25°) ±20° 
-25°) ±20° 
-25°) ± 20° 



Note. — The values for A* indicate posterior mean and stan- 
dard deviation. The ecliptic poles are located at ±(96°, 30°). 



In this paper, we c onsider the seven-yea r WMAP tem- 
perature sky maps (jJarosik et al.l I2010T ) . and analyze 
the V-, W and Q (61, 94 and 41 GHZ, respectively). 
The V- and W-bands are believed to be the cleanest 
WMAP bands in terms of residual foregrounds. Wc 
adopt the template-corrected, foreground reduced maps 
recommended by the WMAP team for cosmological ana l- 
ysis, and impose the KQ85 mask (jGold et al.l I2009D . 
Point source cuts are imposed in all cases. We analyze 
the data frequency-by-frequency, and consider the com- 
binations V1+V2, Q1+Q2 and Wl through W4. The 
noise RMS patterns and beam profiles are taken into ac- 
count for each DA individually. The noise is assumed 
uncorrelated. For d etails on joint Gibbs a nalysis of multi- 
frequency data, see lEriksen et al.l (|2004bD . All data used 
in this analysis are available from LAMBDA. 

The angular resolutions of the V-, W- and Q bands are 
0.35°, 0.22° and 0.53°, respectively. The sky maps are 
pixelized at a HEALPix resolution of iV s id e = 512 with 7' 
pixels. We adopt a harmonic space cutoff of £ m ax = 800 
for the two data sets, probing partly into the noise domi- 
nated regime. However, we do not consider multipoles at 
£ > 400 for the anisotropic part of the signal covariancc 
matrix, in order to minimize the chance of systematic 
effects such as residual point source contributions, beam 
uncertainties or noise mis-estimation to affect our results. 

5.1. Results 

We present the marginal posteriors for the dipolc 
model obtained from the seven-year WMAP tempera- 
ture sky maps, as computed with the method described 
in Section[3J First, in the top row of Figurc|H]we show the 
preferred direction posteriors, P(n|d) for the W- and V 
band data. The joint posterior P(A*,n * |d) is depicted 
in the middle row, while the bottom row displays the 
marginal posteriors,P(„4*|d), and P(n*|d). The results 
are listed in Table [TJ 

The direction of the dipole amplitude is located at 
about (l,b) = (225°, -20°) ± (20°, 20°) in all the 
bands. This corresponds to the direction of the pre- 
viously discovered h emispherical power asymmetry by 
IHoftuft et all (pOOl ) at (l,b) = (225°, -27°), and sug- 
gests that the signal has a cosmological origin. The tilt 
parameter n* is found consistently around 0, implying 
that the dipolc effects are mostly concentrated on large 
scales for £ < 50. The coupling strength of the dipolc 
contribution A* for various bands are listed in Table [TJ 
Note the value of A* is close to zero, there are several con- 
taminating sources that contribute to the posterior. One 
source is the fact that there are degeneracies between A* 



and n* for low values of A*, so the marginal posterior 
have noise- related contributions close to 0. In addition, 
when A* is close to zero, there is a contribution from the 
symmetrical posterior from negative A* that spill over 
to positive values. When keeping this in mind, it should 
altogether be clear that when inspecting the posteriors 
that the amplitude parameter A* is nonzero at about a 
2.0cr significance in the W-, V- and Q bands. 

Since the canonical form of non-commutativity, given 
by Eq.([T]) with 8 a constant in the comoving frame, is 
of specific interest, we have tested this case separately. 
Basically this amounts to fixing the spectral index to the 
theoretical prediction n* = 2 and varying only three pa- 
rameters. In that case we that the data is consistent 
with a vanishing anisotropic contribution and bounds 
A < 0.05. This translates into a coupl e of orders of 
magni tude looser bounds than obtained bv lAkford et al.l 
(2009) fr om the power spectrum alone, since the effect 
is at small scales which we cut out from the analysis 
£ > 400. 

6. CONCLUSION 

In this paper, wc have developed a numerical method 
for investigating traces of non-commutative geometry de- 
rived from field-theoretical implementations of micro- 
causality violation in the seven-year WMAP data. The 
deformation of Lorentz symmetry relevant at inflation 
induces parity-violating modulations of the primordial 
power spectrum, which give rise to a dipole-modulation 
effect in the CMB. The dipolc modulation has cer- 
tain similarities to the previously detected hemispherical 
power asymmetry, and not surprisingly, we reproduced 
the direction of the hemispherical power asymmetry at 
(l,b) = (-225°, -25°) ± (20°, 20°) when analyzing the 
combined seven-year data sets. In addition, both the 
direction and amplitude are stable and nonzero at a la 
level in the W-, V- and Q bands. The tilt parameter n» is 
firmly located around zero, indicating that the seven-year 
WMAP data prefers the dipole modulations to occur on 
large scales £ < 50. While these results are intriguing, 
the significance is still too low to be considered a clear 
detection. This could be due to the fact our parameter- 
ization of the leading order effect may not capture the 
underlying physics to the full effect. In addition, one 
should take into account the higher order multipoles and 
vary also the cosmological parameters. 
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Fig. 6. — The main results from the seven-year W (left) and V (right) band WMAP data. Note that the directions are stable in all bands 
(including the Q-band), and that the anisotropy amplitude At is nonzero at 2<r in all cases. The direction of the previously described 
hemispherical power asymmetry is marked with a red circle. 
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